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We report new results regarding a source of polarization entangled photon-pairs created by the pro- 
cess of spontaneous parametric downconversion in two orthogonally oriented, periodically poled, bulk 
KTiOPCU crystals (PPKTP). The source emits light colinearly at the non-degenerate wavelengths 
of 810 nm and 1550 nm, and is optimized for single-mode optical fiber collection and long-distance 
quantum communication. The configuration favors long crystals, which promote a high photon-pair 
production rate at a narrow bandwidth, together with a high pair-probability in fibers. The quality 
of entanglement is limited by chromatic dispersion, which we analyze by determining the output 
state. We find that such a decoherence effect is strongly material dependent, providing for long 
crystals an upper bound on the visibility of the coincidence fringes of 41% for KTiOPCU, and zero for 
LiNbCX-j. The best obtained raw visibility, when canceling decoherence with an extra piece of crystal, 
was 91 ± 0.2%, including background counts. We confirm by a violation of the CHSH-inequality 
(S — 2.679 ± 0.004 at 55 s~ 1//2 standard deviations) and by complete quantum state tomography 
that the fibers carry high-quality entangled pairs at a maximum rate of 55 x 10 3 s^THz^mW -1 . 

PACS numbers: 03.67.Mn, 03.67.Hk, 42.50.Dv, 42.65.Lm 



I. INTRODUCTION 

A nonlinear medium exposed to an optical field will 
occasionally emit several other photons. The pheno- 
menon is known as spontaneous parametric downconver- 
sion (SPDC), and is frequently utilized for the production 
of photon-pairs. Such a pair can also become entangled 
in a certain degree of freedom if indistinguishability is en- 
sured in all the remaining degrees of freedom. Many suc- 
cessful examples of direct creation of entangled photon- 
pairs Q, 0, , post-selected entangled pairs [j, lail ; an d 
in-fiber generated pairs 0, H, can be given, already 
serving as an indispensable tool for quantum communi- 
cation. 

The source reported here uses two orthogonally ori- 
ented crystals, each emitting pairs of photons of a differ- 
ent polarization than the other. The different pairs are 
made indistinguishable, in our case by single-mode fibers, 
and therefore the individual photons of a single pair be- 
come directly entangled in polarization — an idea orig- 
inally proposed by Hardy [l(| and realized in modified 
form by Kwiat et al. |ll| . One problem with the original 
realization is that the crystals cannot be made too long, 
since the non-colinearity makes the two emission-cones 
non-overlapping. Another problem is that the crystals 
used generally emit into many spatial modes, which is 
not suitable for fiber- coupling. U sing pe riodically poled 
crystals via quasi-phase matching |l2lll3LIT^ |. it has been 



'Corresponding author. Electronic address: daniellj@kth.se 
URL: http://www.quantum.se 

t Current address: Department of Physics, University of Toronto, 
Toronto M5S 1A7, Canada 

•f Current address: Department of Physics, University of Chicago, 
Chicago, IL 60637, USA 



shown that colinear emission can be achieved very close 
to a single mode }l5j | (even in non-waveguiding struc- 
tures), providing much greater overlap in the emission. 
Such a configuration also allows non-degenerate wave- 
lengths to be generated. 

Some desirable properties of sources to be used for 
quantum communication include: i) a high probability 
of photon-pairs to be collected into optical fibers; ii) a 
minimum number of false coincidences; iii) wavelength 
combinations that either suit efficient detection, match 
atomic transitions, or are well transmitted over long dis- 
tances; iv) a narrow bandwidth that limits the effects 
of fiber dispersion (~ GHz) or can address atoms 
(~ MHz) ; v) a long coherence length that limits the need 
for precise interferomctry; vi) small jitter in arrival-time 
of photons; vii) perfect correlations in all bases; and, ide- 
ally, viii) the source being compact enough to be put in 
a box, carried out of the lab, and be used, e.g., for quan- 
tum key distribution (QKD). Furthermore, for maximum 
security in QKD a strong requirement is to have neither 
more nor less than a single photon per gate pulse. In this 
respect, photon-pair sources have been shown to be good 
candidates compared to weak coherent pulses, potentially 
fulfilling properties i) and ii). Equally imperative for se- 
curity in Ekert's scheme |17| is property vii), which ex- 
presses the wish for high visibility of entanglement in the 
presence of background detection, which implies the need 
to minimize dark counts and false coincidence counts. 

In this work, we extend our previous results |l8j re- 
garding a PPKTP-based two-crystal source and try to 
address some of the anticipated features above. By emit- 
ting at non-degenerate wavelengths, the source exploits 
the highly efficient Si-based single-photon counters avail- 
able in the near-infrared region and the low attenuation 
in fibers at telecom wavelengths. The shorter wavelength 
also matches the transmission bands of alkaline atoms, 
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FIG. 1: The source consists of two PPKTP crystals placed one 
after the other; the first creates a vertically (V) and the second 
a horizontally (H) polarized field. The laboratory coordinate 
system is drawn, as well as the crystal axes, X, Y, and Z, 
which refer to the polarization of the incoming and outgoing 
electromagnetic fields. 



which makes the source suitable as part of a quantum 
memory [l9L |20|. For our crystal configuration, we show 
how effects like chromatic dispersion enter the picture 
as problems to be dealt with. The source has been 
optimized for coupling into single-mode fibers following 
Ref. |l5l| . where one can also find motivations for us- 
ing long crystals to achieve a narrow bandwidth. An 
early example of a non-degenerate source is Ref. [2lJ , us- 
ing energy-time entanglement. One reason for utilizing 
energy-time entanglement is to overcome the strong de- 
coherence mechanism of polarization-states over fibers, 
and, for the same reason, we propose a scheme that com- 
bines time-multiplexed encoding on the telecom wave- 
length side |22j | with polarization on the near-infrared 
side, altogether realizing a sort of hybrid-coded entan- 
glement. 

The article is organized as follows. In section II, we 
describe the main characteristics of the source. In sec- 
tion III, we derive the quantum state emitted by the two 
crystals in terms of frequency and polarization degrees 
of freedom, based on the quantum state of a single crys- 
tal derived in the Appendix. Following that, in section 
IV, we briefly show how to compensate for the effect of 
chromatic dispersion in the crystals, so as to assure indis- 
tinguishability, and, in section V, we present our experi- 
mental results showing the quality of the source, includ- 
ing results on quantum state tomography. In section VI, 
we discuss the future directions of a hybrid-coded source, 
and we end with a summary in section VII. 



II. A SOURCE OF POLARIZATION 
ENTANGLEMENT 

The source is depicted in Fig. ^ an d consists of two 
orthogonally aligned bulk crystals placed one after the 
other. They each have the dimensions 3 x 4.5 x 1 mm 
(X, Y, Z) , of which the second dimension defines the 
length, L = 4.5 mm. The crystals are made of potassium 



titanyl phosphate, KTiOPC>4, and are periodical poled 
with the period A = 9.6 /^m, chosen such that we have 
phase-matching for the signal at a wavelength of 810 nm, 
and the idler at 1550 nm, for a temperature T = 111°C 
determined by the Sellmeier equations of KTP (2^, |3 . 
The crystals are pumped by monochromatic and con- 
tinuous wave laser light (p) at a wavelength of 532 nm, 
which is propagating in a Gaussian TEMoo mode along 
the z-axis, producing a signal (s) and idler (i) field in the 
same direction and with the same polarization as the Z- 
component of the pump field (Z V Z S Z{). Fig.^defines the 
laboratory axes and the crystals' optical axes X, Y, and 
Z, oriented as shown. Both crystals will generate down- 
converted light if the pump polarization is oriented at 
45° to the horizontal (H) and vertical (V) axes. Follow- 
ing |l "J| , we have optimized the focusing of the pump and 
the fiber-matched modes using the parameter £ = L/zr, 
where L is the length of the crystal and zr is the Rayleigh 
range, such that the maximum amount of emission that 
is generated is collectible into single-mode fibers. The op- 
timal values for our configuration arc £ p = 1.3, £ s = 2.0, 
and = 2.3, respectively. 

The use of single-mode fibers to collect the light will 
erase all spatial information that reveals from which crys- 
tal the photons came, except for the polarization degree 
of freedom. Therefore, each of the beams will interfere 
in the diagonal basis and get entangled in polarization. 
(Note that the spatial information is partly correlated 
with frequency via the phase-matching condition, and 
that indistinguishability could also be achieved via fre- 
quency filtering.) The resulting state is the Bell-state, 

|^) = ^(|V} s |V) i + e^|H) s |H) i ), (1) 

with a relative phase if that we can control. As in most 
cases, it is required that the probability of creating more 
than a single pair within a time determined by the co- 
herence time of the photons, or the detector gate-time, 
whichever is longer, is negligible, and for moderate pump- 
powers and relatively short gate-times or wide band- 
widths, this probability is very small, but not vanishing. 
Assuming a Poissonian distribution, the probability be- 
comes P n >2 = 1 — (l + m)e _m , where m = At g /3P p Xp/hc 
is the mean photon number in a single random gating. 
For a typical detector gate-time At g = 5 ns, pump-power 
P p = 540 //W, and conversion efficiency (3 = 3 x 10~ 10 
we get m = 2 x 10~ 3 and P n > 2 = 2 x 10~ 6 . 

Fig. [21 will serve as an illustration of the problem of 
optimizing the focus of the pump-mode, and the fiber- 
matched modes with respect to two crystals. As a com- 
promise, the pump-beam is focused at the interface be- 
tween the crystals, in the anticipation that the profile of 
the generated emission exactly trails the profile of the 
pump-beam. However, numerical simulations with the 
software developed in 0] show that the waist of the 
emission will be shifted towards the center of each crystal, 
so that neither the vertically nor the horizontally polar- 
ized photons will couple perfectly into the fiber simulta- 
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FIG. 2: The figure shows a Venn diagram, which is used 
to illustrate the single coupling efficiencies 7 S and 71, pair 
coupling 7c, and conditional coincidences fi B ^ and jUii 8 as a 
fractional number representing the area of a set. The elements 
contained in a specific set represent photon pairs that are 
coupled into a fiber taken from the universal set of pairs, 
fi p , which contains all pairs generated from the crystals (V 
or H) within the bandwidth of the detector filter AA. The 
set fi p is normalized to unity and represents perfect coupling. 
Maximum overlap of all sets is needed to generate the best 
entanglement in the fiber, which is represented by the darkest 
shaded area in the diagram (the union of all sets). 



neously. The figure shows the different types of coupling 
efficiencies represented as sets in a Venn-diagram, where 
each element of a set represents a photon pair generated 
by the crystals in some spatial mode. That is, the collec- 
tion of all elements within each set defines which pairs are 
coupled into the fiber for some specific focusing condition, 
in such a way that the coupling efficiency corresponds to 
the total area of the set. The problem can be described in 
two parts: first, the need to overlap the matching modes 
of the signal and idler, represented by the coupling effi- 
ciencies 7 S and 71, for each polarization separately (i.e. 
by optimizing the pair coupling 7 C = /ii| s 7 s , via the con- 
ditional coincidence (J,i\$), and second, the need to over- 
lap the vertically, 7 V , and the horizontally, 7 H , polar- 
ized photons for both the signal and idler, ft is only 
in the intersection of all sets where entanglement exits, 
and any detection of photons outside of this set will limit 
the visibility in the ±45°-basis (denoted here D/A-basis) 
by contributing to a mixed state. This picture is valid 
for many types of sources, and we believe that the cou- 
pling efficiencies in many cases in the literature are esti- 
mated in an incorrect way, as it is important to note that 
7c 7^ 7s 7i (especially in non-degenerate regimes). By this 
short discussion (see plj for a comprehensive discussion) , 
we hope to have illustrated that it is not necessarily best 
to optimize each arm individually to find the greatest co- 
incidences, but rather, to simultaneously optimize both 
arms. 

As we have mentioned, the different polarizations need 
to interfere, and therefore a major concern is that they 



FIG. 3: Color online. The figure illustrates the effects of chro- 
matic dispersion resulting from the strong non-degeneracy of 
the signal and idler photons, with group refractive indices n g , s 
and n gi i. a) A pair created at the end-facet of the V-crystal 
(black dots) will non-interactively pass through the H-crystal, 
after which the signal and idler wavepackets become separated 
by a time ti = (n^ s — n^^L/c before hitting the detectors. 
As seen from the detectors' viewpoint, for the same pair to 
instead have been created in the H-crystal (black dots), its 
wavepackets would necessarily need to have separated by the 
same amount, given by t[ = (n? s — ngi)b/c, in order to inter- 
fere with (i.e. overlap with, or be indistinguishable from) the 
first case. (The superscripts refers to different polarization- 
axes.) Any pair created within the gray area (between black 
and white dots), are separated by a time r > n = t[, and 
will always find a corresponding position in the orthogonal 
crystal to interfere with, according to the detectors; however, 
all the pairs from the white area of either crystal will be dis- 
tinguishable in time from any pair of the other crystal, i.e. the 
wavepackets are non-overlapping due to different dispersions, 
and contribute therefore to a mixed output state, b) If we 
put a birefringent plate of thickness d in one of the arms, the 
time-separation for a pair created at the end-facet of the V- 
crystal is reduced to r 2 = (n* s — n^ ; )L/c+ (1 — n^ c )d/c < n, 
for which some d equals the time-separation of an interfering 
position at the end-facet of the H-crystal, T2 = — (1 — n° c )d/c. 
Consequently, all H and V-pairs now show "self-interference" , 
and a pure output state is created. Note that two pairs cre- 
ated within the coherence time of the pump (which needs to 
be longer than L) are always coherent, ignoring dispersion. 



are not distinguishable by time information, noting the 
limited extent of the photon wave-packets. For long crys- 
tals, the photon pairs will separate by chromatic disper- 
sion, due to the very different group velocities between 
the strongly non-degenerate signal and idler. This will 
occur to a degree that is different for pairs created in 
the first crystal than for pairs created in the second, be- 
cause the pairs from the first crystal also need to pass 
through the second. The differences in group velocities 
between signal and idler are not the same for light polar- 
ized along the iT-axis and the X-axis, implying that not 
all pairs, created along the length of cither crystal, will 
find any (possibly) generated pairs to interfere with from 
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the other crystal. Some photon-pairs will therefore be 
distinguishable by temporal information. See Fig. for 
an illustration. We would like to point out that, while 
this chromatic two-photon dispersion effect is reminis- 
cent of the "two-photon dispersion" effects discussed in 
PJ or pfij . it does not have the same origin, although 
the current effect can also be compensated for by an ex- 
tra piece of crystal. The chromatic effect comes as an 
disadvantage when placing the crystals adjacent to each 
other, and could inprinciple be avoided by an "interfer- 
ometric" solution |y, [l2| , in which the pump beam splits 
into two separate arms, impinges onto each of the crys- 
tals, or onto a single crystal but in opposite directions, 
and recombincs on a beam-splitter. Still, we believe the 
current solution requires fewer optics, is easier to align, 
and can be made more compact. 

The previous discussion gave a limited, although intu- 
itive, understanding of the origin of a mixed state, but, as 
we will show in the next section, a mathematical deriva- 
tion will give additional insights into how the effect of 
decoherence is affected by the group velocities. 



III. THE TWO-CRYSTAL TWO-PHOTON 
QUANTUM STATE 



In this section, we derive the output state from the two- 
crystal source in terms of the frequency and polarization 
degrees of freedom \e)®\Xij)> where i, j = {1 = "V", 2 = 
"H"} denotes the polarizations. Emission from each of 

sentcd by 



and IX22), respectively, according to Eq. (| A. 1 5|) of the 
Appendix and Fig.^ As just described, the vertical light 
will be subject to dispersion upon its passing through the 
second crystal. We will formulate this mathematically by 
introducing a unitary transform acting on the states. The 
cigencquation which describes the transformation Ul on 
the state of the first crystal, when it passes through the 
second crystal, is 



will give us the combined two-crystal two-photon state, 



U L \xxi) = e" v " BiJ 1 '""•'IXii 



i(k B L+hL) j 
i(nf u B+n^uioi+(nf tB -nf >i )e)L/c\^\^ 



where the length of the crystal, L, enters the phase term, 
together with the frequency e. With reference to the Ap- 



pendix B and Eq. 
state of each crystal as 



we can then express the output 



l*n) = ^ JdeU L U(e)\e)®\ X n), 
|*22> = 4 fde U(e) e m >^ c |e> ® | X22 >, (3) 
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where an extra phase-term has been added to the pump 
field in the second crystal due to the pump field passing 
through the first crystal, U(e) is the state amplitude, and 
B is a normalization constant. The sum of these two kets 
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+ 1*22 

1 
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de 



U L U(e)\e)®\ X 



+ U(e)e m >- L / C \e) ® | X22 } 



I f 2 

- Jde CijVijie) |e) <g> \xij), 



(4) 



where we have introduced Vh(e) = -^UlU(€), V^e) 
±U(e)e m p u> * L/c , and the coefficients c v = I/V2 for i 
j, and Cij — for i 7^ j, normalized such that |cn| 2 
|c 22 | 2 -f. 

We can now form the two-photon density matrix 
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\i JJ dedi E Vy(e)tS(c) |e)(e| 



B 



i,j,k,l—l 



\Xij)(Xkl\, 
(5) 



from which we would like to remove the frequency in- 
formation. For that, we need to note that we could, in 
principle, measure the frequency of the photons at a res- 
olution much smaller than the bandwidths of the filters. 
The resolution is given by a wavelength bandwidth AA res , 
which is set by the timing information Ai ga t of the de- 
tectors (AA res = A 2 /cAi gato < 8 pm for At gatc > I ns). 
Therefore, it is appropriate to take the partial trace over 
the frequency mode: 

00 

p = Tr e [| e ')< e '|p e ]= Jde'(e'\p'\e') 

— OO 

= E ^c* kl Jde'V l3 {e')V^e')\ Xl3 ){XM\- (6) 

i,j,k,l— 1 

Let Pijki denote the elements of the density matrix, of 
which the only non-zero ones become 

Pii22 = ciic^ Jde'U L U(e') U* {e')e~ m >- L ' c 



2 h 2 B 2 

x [de>\A a (e>)\*\A(e')\ 2 



X e p p ' e 

'Le' 



x sine 



P2211 



-^«s-<i) 



and 
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FIG. 4: Color online. The off-diagonal term, Eq. of the 
generated density matrix plotted versus crystal length, which 
corresponds to the visibility of entanglement via the relation 
V = 2/91122- Solid line: PPKTP-crystal with 10 nm idler 
filter. Dots over solid line: PPKTP with 1 nm idler filter. 
Dash-dotted line: PPLN or PPMgOLN with 10 nm idler filter. 
Dashed line: PPKTP at optimal fiber coupling using the idler 
fiber's own filtering. Diamond: experimental value for L — 
4.5 mm using PPKTP and with a 10 nm idler filter. Solid 
point: experimental value for L — 4.5 mm using PPKTP and 
with a dispersion canceling calcite plate of thickness d = 0.86 
mm. 



The off-diagonal element, which describes the degree of 
coherence in the entangled state, can be further simplified 
and identified as a Fourier transform: 

p 1122 — - e ~ in i "p L / c e K n f u °a+ n ?<*>oi)L/c 

x fde' g(e')e^' sinc 2 (^e'), (9) 



where 



g(0 = ^4§^ \A s (e')\ 2 \Me')\ 2 , 



h 2 B 2 
z z 



(n 



;,s - n ti) L / C - 



(10) 



In Fig. QJ we have plotted the result of Eq. versus 
the length of the crystals using different crystal materi- 
als to generate 810 and 1550 nm. We observe that the 
dispersion in long, periodically poled LiNbC>3 (PPLN) 
crystal materials completely suppresses the P1122 term, 
and thereby the entanglement. For PPKTP this is not 
the case, and if we search for P1122 in the limit of an 
infinitely long crystal we find that 



lim I/0H22 1 




(11) 



which, for PPKTP leads to p n22 = 0.203, implying still 
a visibility of entanglement (i.e. of the second-order in- 
terference fringes) of 40.6%. The different results stem 



from the material-specific relation between t x and r z . If 
the material is more strongly dispersive for polarizations 
along the Z-axis than the X-axis, then the off-diagonal 
term will be bounded below by a non-vanishing value; 
otherwise, the off-diagonal term will approach zero. As 
expected, we found that all numbers increase as we go 
closer to having degenerate wavelength pairs. We also 
note that the bandwidth of the frequency filter affects the 
shape of the curve; a narrower bandwidth increases the 
extent of the temporal coherence and provides a greater 
overlap between wave-packets, leading to an arbitrarily 
increased pn22- Emission that is optimally coupled into 
single-mode fibers will automatically be filtered also in 
frequency, since the frequency is correlated to spatial in- 
formation via the phase-matching conditions [l5j , and for 
long PPKTP crystals, in such a case, the minimum value 
of P1122 equals 0.266 (V = 53.2%). 



IV. DECOHERENCE CANCELLATION 

We will now briefly show how the pure state, in 
Eq. Q), can be fully regained, for generation in long crys- 
tals, by inserting a highly birefringent crystal plate into 
one of the arms. The eigenequations for each polarization 
state propagating through such a crystal plate become 



Uc\x 



U C \X2j) 



e i«iJoi-nJ iC e)d/c|^^^ 

e \X2j) 

e i(n^uj 0i -nl a e)d/c^^ 



(12) 



with the density matrix after the plate becoming 

p^d) = U c p f ul. (13) 

Magnitude off-diagonal term, p 
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FIG. 5: Color online. The off-diagonal term, Eq. H14II . plot- 
ted versus the crystal lengths, L, and the thickness, d, of a 
dispersion-canceling calcite plate, using PPKTP with a 10 nm 
idler filter. The dashed line represents perfect cancellation. 
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Repeating Eqns. J5J to {UJ, we arrive at 



P1122 = - e~ in p U ? L / C e i( - n * <"o s +n*w i)£/ c 

v i(ri°w i-«c w Oi)d/c 



x j&e' g{e')e t{T x- K)€ ' sinc 2 (^e'), (14) 

where g(e'), t x , and is defined by Eq. I|1U|) . and where 

k = (n° c - nl iC )d/c. (15) 

Now, if d is chosen such that k = t x , it means that we 
have perfectly canceled the decoherence and retrieved a 
pure state. Hence, 



1 



P1122 — P2211 ^ P1111 — P2222 



(16) 



Note that, by adjusting d and tilting the plate (affecting 
tp) our source can prepare any arbitrary mixed state of 
the kind p = V\^ v ){^ v \ + (1 - V)p m , where p m = 
^(IXiiXxnl + IX22XX22I), and V is the visibility. Fig.0 
shows a plot of P1122 versus L and d. 



V. EXPERIMENTAL RESULTS 

The experimental setup used when characterizing the 
source's output state is shown in Fig. |SJ As a pump, we 
use a frequency-doubled Nd:YAG laser emitting approxi- 
mately 60 mW in the TEMoo mode at 532 nm, which can 
be variably attenuated. Its M 2 value was measured to 
be 1.06. After a band-pass filter (BP532) that removes 
any remaining infrared light, we "clean up" the polariza- 
tion using a polarizing beam-splitter (PBS). The polar- 
ization is controlled by a half-wave plate (HWP) and a 
quarter-wave plate (QWP) in front of the crystal. The 
pump beam is focused onto the crystal using an achro- 
matic doublet lens (/ p = 50 mm), which introduces a 
minimal amount of aberrations, so as not to destroy the 
low M 2 -value. The QWP is set to undo any polarization 
clliptisation effects caused by the lens, and fluorescence 
caused by the same lens is removed by a Schott-KG5 
filter (SP). 

The next components are the two PPKTP crystals, 
which are heated in an oven to a temperature T w 100°. 
After the crystals, we block the pump light with a 532 nm 
band-stop filter, and the signal and idler emission is fo- 
cused by achromatic doublet lenses. To separate the 810 
nm and 1550 nm emission, we use a dichroic mirror made 
for a 45° angle of incidence. The first lens (/ s ; = 30 mm) 
is common to both signal and idler, and its task is to refo- 
cus the beams somewhere near the dichroic mirror. The 
next two lenses (/ s = 60 mm and /i = 40 mm) collimate 
each beam, which are then focused into the fiber tips 
(with the mode field diameters being MFDgio = 5.5 /im 
and MFD1550 = 10.4 /im) using asphcrical lenses with 
/ = 11 mm. Next, we use quarter-wave plates (QWP), 



half-wave plates (HWP), and polarizing beam-splitters 
(PBS) in each arm to analyze the state. In the idler 
arm, we also place the tiltable cancellation plate, which 
is made of calcite. In front of the fiber couplers, we have 
first Schott-RG715/RG1000 filters to block any remain- 
ing pump light, and then interference filters (BP) of 2 
nm and 10 nm bandwidth at the 810 nm and 1550 nm 
side respectively. The detectors used are a Si-based APD 
(PerkinElmer SPCM-AQR-14) for 810 nm with a quan- 
tum efficiency ?7 S = 60% and a homemade InGaAs-APD 
(Epitaxx) module for 1550 nm with rji = 18%, gated with 

5 ns pulses. To avoid afterpulsing effects, the InGaAs- 
APD is used together with a hold-off circuit (10 /zs) for 
all of the measurements. The pulses were generated us- 
ing a digital delay generator (DG535) from SRS, with a 
maximal repetition rate of 1 MHz, and a trigger dead 
time of 1 /xs. 

We have used a spectrograph (SpectraPro 500i, ARC) 
to measure the bandwidth of the signal emission using 
a single-mode fiber without any filter; see Fig. [7| The 
bandwidth was found to be 4 nm for the V-crystal and 

6 nm for the H-crystal. The results in suggests that 
the effective lengths of the crystals being poled must then 
be 3 mm and 2 mm, respectively, but also that the 2 mm 
crystal should give w 55% of the photon-rate of the 3 
mm one. Experimental agreement is good, as we saw 
the H-crystal giving half the rate of the V-crystal. When 
measuring, we refocused the fiber coupling for each crys- 
tal to find maximum counts, while keeping the pump po- 
larization exactly at 45°. As described in connection to 
Fig. [21 the best tradeoff when collecting from both crys- 
tals simultaneously is to set the focus of the pump mode 
and the fiber-matched modes at the intersecting faces. 
Experimentally, however, in order to produce as pure a 
Bell-state as possible, we needed to balance the rate of 
each crystal, which we did by shifting the fiber-matched 
focus a bit closer to the H-crystal and by turning the 
pump-polarization slightly towards H. (The focus point 
was moved by turning the focusing knob on the fiber 
coupler.) In this way we allowed lower coupling efficien- 
cies than the maximum attainable. The focusing condi- 
tions achieved with available lenses were, £ p = 2.1 for 
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FIG. 6: Experimental setup used to measure the density ma- 
trix. PBS: polarizing beam-splitter, HWP: half- wave plate, 
QWP: quarter- wave plate, SP: short-pass filter, BP: band- 
pass filter, SMF: single-mode fiber. 
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TABLE I: Three runs at different alignments and pump powers, showing the coupling efficiencies, photon rates in fibers, 
conversion efficiency, and the production rate of the system. 
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FIG. 7: Color online. Spectrogram of signal emission inside 
single-mode fibers without interference filters. The crosses 
(AA S = 6 nm) and circles (AA S = 4 nm) represent experi- 
mental data for the H and the V crystal, respectively, to be 
compared to theoretical predictions for a 2 nm (solid line) and 
a 3 nm (dashed line) long crystal. Also shown is the downcon- 
version spectrum in free space, experimentally (dots) and in 
theory (dotted line), demonstrating the fiber's own filtering. 



the pump mode, £ s = 3.2 for the signal's fiber-matched 
mode, and £i = 2.5 for the idler's. 

With this configuration, we obtained the results 
showed in Tabled In each column of the table, 7 S repre- 
sents the signal's single coupling efficiency, 7; the idler's, 
7c the pair coupling efficiency, fj,u s the conditional coinci- 
dence, and a the correlation efficiency, which corresponds 
directly to /z;| s but includes a compensation for the 35% 
transmission of the 1550 nm filter and the 85% trans- 
mission of the 810 nm filter. The singles photon-rate in 
the signal fiber, R s , and the idler i?;, were both derived 
from detected raw counts. The total generated rate R p of 
pairs before fiber coupling was estimated from detected 
counts using a multimode fiber. The pair rate in the 
fibers, i? c , was deduced from the above efficiencies and 
the detected raw coincidence rate, with accidental counts 
subtracted by assuming that i?j originates from a Poisso- 
nian distribution at random gating |25j |. The conversion 
efficiency f3 is the fraction of pump photons converted 
into signal and idler pairs, leading to a pair production 
rate R% Iod , which equals 5 x 10 3 s -1 THz- 1 mW _1 at the 
pump power P p = 60 mW and with the idler detector 
gated at 585 kHz. (The production rate is the pair rate 
normalized to the wavelength bandwidth in THz and the 



pump power in mW.) The second row of Tabled shows 
similar results for a lower pump power and an idler gate 
rate of 91 kHz. We also took measurements without any 
interference filter at the idler side (but still with a 2 nm 
filter at the signal), with the results shown in the third 
row of Table U for P p = 540 fiW and with a gate-rate 
of 57 kHz. The results are improved, not because of the 
lower power, but because of a simultaneous optimization 
of the arms in order to maximize j c . The table shows how 
7i decreases in the process. The correlation efficiency a 
now includes the estimated transmission-loss of the optics 
at the idler side, and a correction factor for the unequal 
filtering between signal and idler (the idler fiber itself pro- 
vides a frequency filtering of AAi = 14.7 nm). The best 
conditional coincidence is [i^ s = 0.27, and the conver- 
sion efficiency, (3 = 3 x 10~ 10 , was possibly improved by 
aligning to a more homogeneously poled area of the crys- 
tals. We believe that the pair production rate, 55 x 10 3 
s _1 THz _1 mW _1 , is one of the highest yet reported for 
polarization entangled photon pairs generated in crystals 
and launched into single-mode fibers. Frequency filtering 
at a narrow bandwidth of 50 GHz would imply 3 x 10 3 
s^mW -1 of pairs in the fibers. Besides, for narrow fil- 
tering, the photon flux has been shown [l5| to be oc L\/L, 
and so, by using longer crystals (L = 50 mm) we could 
still reach 20 s _1 mW _1 at a 10 MHz bandwidth, which 
is the bandwidth regime of e.g. Rb-atom based quantum 
memories. As a comparison, we have derived numbers 
using data available for some other experiments, among 
which the best include Fiorentino et al. [l^ 1 wn0 seem to 
have 22 x 10 3 s _1 THz _1 mW _1 of pairs being generated 
by two 10 mm long crystals into free-space; Konig et al. 
|20j . who claim to have 300 x 10 3 s _1 THz- 1 mW _1 pairs 
from two 20 mm long crystals into fibers; and Li et al. 
0, who seem to have an exceptional value of 4.3 x 10 6 
s _1 THz _1 mW _1 pairs generated directly inside a non- 
linear fiber. 

Fig. [5] shows the visibility curves obtained, with and 
without subtraction of background counts, including false 
coincidences. Note that the number of "accidental" co- 
incidences increases with the pump power, as the prob- 
ability of more than a single pair to arrive within the 
gate-time of the detector increases, as shown in section 
II. 

We have also measured a violation of the CHSH- 
inequality |27j by taking measurements of the 
coincidence-rate functions 
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FIG. 8: Color online. The plot shows the raw visibility 
curves obtained at a pump power P p = 4.5 mW, gate-rate 
91 kHz, without subtraction of background counts (includ- 
ing accidentals) shown at the bottom. Each curve corre- 
spond to a different polarization setting of the signal arm, 
(j>Y = -37r/16, <j>f = 7r/16, <ff = -it/16, and <j>£ = 3tt/16, 
as indicated the inset. Fitting of the collected data gives 
V H = 99.6 ± 0.2%, VV = 99.2 ± 0.2%, F D = 92.7 ± 0.2%, 
and Va = 94.2 ± 0.2%. When background counts were not 
subtracted, we obtained the visibilities Vh = 95.6 ± 0.2%, 
Vv = 96.2 ± 0.2%, Vd = 89.6 ± 0.2%, and Va = 90.9 ± 0.2%. 



where i, j = ±1 denotes the four combinations of mea- 
surable output-arms of the two PBS:s in the signal and 
idler, Vi j is the corresponding visibility, and cf) s , fa u are 



'hi 

the angles of the HWPs. 
comes 



The correlation function be- 
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Ri,i — Ri,-i — R-i,i + R-1,-1 
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V"cos(40 s + 4 



R-i, 



R- 



(18) 



where V = Vi t x, assuming fully equal rate functions, so 
that we can rely on measurements taken at only one of 
the output arms. Entanglement is present iff the CHSH- 
inequality is violated, 

S = E(ti,$) + Etflti) + \E{<Plti) - E{<Pl, fat)\ < 2, 

(19) 

where the correlation function is to be measured at the 
following pair of angles: fa] — — n/16,fa{ = and 
fal = n/16,(j>f = 7r/8. The parameter S can reach the 
maximum value of 2\/2, corresponding to 100% visibil- 
ity, and it is well known that the average visibility needs 
to be > 71% to violate the inequality, if the state is 
subject to equal decoherence in all bases. In our case, 
the state decoheres in the H/V-basis, while maintaining 
nearly perfect visibility for the H and V settings. Let fa] 
represent the H/V-basis, and fal the D/A-basis. Further- 
more, let Vh.v an d Vd.a represent the visibilities in each 



respective basis. We get 

S = V* H ,vcos(-| + 0) + Vk,vcos(-| + |)+ 

7T 7T 7T 

Vb,A cos (j + 0) - Vt>,a cos ( 4 + 2 < 



V2(Vk.v + Vd, 



(20) 



which shows that, for Vh,v = 100%, the requirement is 
Vd,a > 41% for a violation of Eq. {TjjJ. 

A direct measurement of S at the above angles yields 
S = 2.679 ± 0.004 at a pump power P p = 60 mW, after 
subtraction of accidental counts (gate-rate was 585 kHz). 
The CHSH-inequality was violated by 177 standard de- 
viations in 10 s, or 56cr s s" 1 / 2 . To our knowledge, this is 
one of the highest reported to date; only Kurtsiefcr et al. 
exceeds this rate, with 148(7 S s -1 / 2 . Other examples 
of good results can be found in 

[HI (50ct s s" 1 / 2 ) and 
in 01 (38c s s -1 / 2 ). We have re-derived these numbers 
using available data, in the hopes of having created di- 
rectly comparable normalized numbers. The derivation 
was made as follows. Assuming no fluctuation of the rate 
other than that originating from Poissonian-distributcd 
single-photon detections, the standard deviation of the 
coincidence rate Rij becomes a R = \J i? max /2/ ^/Tr, 
where i? max is the peak coincidence rate, a/ R max /2 is the 
standard deviation of the average photon-rate, Tr is the 
integration-time in seconds, and the central limit theo- 
rem is used to sum over time. According to Eq. I|18(l . 
the standard deviation of the correlation function be- 

1X , and by Eq. ltT§|) we have 
eucTr, such that S = S m ± cr^ 



comes o E = \fi<r R /2R 



>E " 2/ V^-RmaxJ-Ki ouiai uiiau u — u m _i_ tig, 

where S m is the measured value over Tr seconds. Thus, 
the normalized "speed of CHSH violation" becomes 



S m — 2 (S m — 2)^/2R u 



-1/21 



(21) 



which only depends on the maximum rate and the mea- 
sured value of S. If the accidental counts are not sub- 
tracted from the coincidence counts, we instead mea- 
sure the value S = 2.6283 ± 0.0102 (P p = 4.5 mW), 
with the CHSH-inequality being violated by 19a s s -1 ^ 2 , 
showing that we truly have a high degree of entangle- 
ment launched into the fibers. This is important in 
entanglement-based quantum key distribution (QKD) 
systems that do not allow a subtraction of the back- 
ground. Rather, any accidentals will increase the quan- 
tum bit error rate (QBER) and reduce the final bit rate, 
equivalently degrading the system performance. 

Following Ref. [2^], we have made a complete tomog- 
raphy of the state, with the resulting density matrix be- 
coming 



Pcxp 



0.5197 -0.0237 0.0300 0.4573 
-0.0237 0.0069 0.0146 -0.0114 
0.0300 0.0146 0.0010 
0.4573 -0.0114 0.0010 0.4734 



+ i 



0.0628 -0.0150 0.0720 
-0.0628 -0.1107 0.0206 
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FIG. 9: Color online. Experimentally determined density 
matrix, p C x P , (real and imaginary parts) obtained by quan- 
tum state tomography on the generated polarization entan- 
gled state (1 = "V" and 2 = "H"). (Pump power, P p = 4.5 
mW). 



which is also plotted in Fig. El Recall that the off- 
diagonal element, P1122 = 0.457, corresponds approxi- 
mately to the visibility in the D/A-basis, V « 2pn22 = 
0.915, which is indeed close to the measured visibilities. 
When applying the density matrix to Wootters's entan- 
glement of formation measure p9l | , we get the value E = 
0.56. The entanglement of formation equals unity for a 
pure Bell-state, as do the fidelity, F = ($ ¥, |pexp|*' p ), 
which is found to be 0.95 for the generated state. 



VI. FUTURE DIRECTIONS: A 
HYBRID-CODED ENTANGLEMENT SOURCE 

In order to motivate the usefulness of the source, we 
provide in Fig. I1UI a complete setup for quantum com- 
munication (e.g. QKD). The scheme, which is under im- 
plementation, uses long crystals (2 x 50 mm) in order to 
achieve a bandwidth of < 80 GHz, which means higher 
production rates and less dispersion in combination with 
a telecom Bragg grating as dispersion compensator. For 
long crystals, the optimal focusing is weaker, which leads 
to a more compact source with fewer collimating lenses 
placed at closer distances to each other. Furthermore, 
improvement of the conditional coincidences as well as 
the size of the source can be achieved by minimizing the 
number of components, each of which contribute to loss. 

In Bob's arm, the polarization information is converted 
into time information in order to avoid the polarization 
dispersion in standard telecom fibers. (For a thorough 
review on photonic qubits, please refer to H3-) A polar- 
izing beam splitter sits in an unbalanced Mach-Zchndcr 
interferometer, directing vertical photons into the long 
arm and horizontal into the short. The vertical pho- 
tons are rotated to horizontal before the photons in both 




FIG. 10: Scheme to create hybrid-coded entanglement. BS: 
beam splitter, PBS: polarizing beam splitter, HWP: half-wave 
plate, BP: band-pass filter, BG: Bragg grating, FM: Faraday 
mirror. 



arms are recombined on a fiber-based beam-splitter and 
sent to a Bragg grating. The result is a time encoded 
qubit with all polarization information erased. (To our 
knowledge, there is no way to erase this information pas- 
sively without having to accept 50% losses in the un- 
used arm of the beam-splitter, which is an disadvantage, 
but could also be turned to an advantage by introduc- 
ing a third party Charlie.) The resulting state becomes 
|$') = l/V2(|V) s |L)i + e^|H) s |S)i), where L denotes 
the long arm and S the short arm. On Bob's analyzer 
side, there is an unbalanced all-fiber Michelson interfer- 
ometer with a single beam-splitter to decode the qubits. 
The interferometer uses Faraday mirrors, which reflect 
the light in such a way that the polarization is exactly 
orthogonal when the photons arrive a second time at the 
beam-splitter to interfere, and thereby avoids the need 
for polarization controllers [3lT |. The phase information 
of the qubit defines a complementary basis to time, and 
for that information to remain, the path length differ- 
ence between the short and the long arm needs to be 
exactly matched to that of the preparing interferometer, 
requiring both interferometers to be temperature stabi- 
lized. However, longer coherence length of the emitted 
photons (an effect of narrow bandwidth) will effectively 
relax these requirements. One advantage of the above so- 
lution is that the preparing interferometer has translat- 
able fiber couplers inside the interferometer, which sim- 
plifies their mutual alignment. Also, we avoid the possi- 
bly difficult alignment of three interferometers, as Alice 
adheres to polarization coding. Another important con- 
dition for the qubits to remain coherent is that the delay 
between two consecutive pulses is short enough 5 ns) 
that they experience the same phase shift due to vibra- 
tions and temperature fluctuations when traveling over 
the fiber. On Alice's side, the analyzer realizes a standard 
polarization decoder. Note that the H/V or D/A-basis 
is randomly chosen by the first beam-splitter, just as at 
Bob's side, which implies that there is no need for any 
active devices. Note also that there exists the possibility 
to delay the outputs of each detector arm on Alice side 
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and combine into different time-slots for detection with 
a single detector, instead of four, which may reduce the 
need for space. 



VII. SUMMARY 

In this article, we have presented work on a two- 
crystal source that uses PPKTP for the production of 
polarization-entangled photon-pairs in a single spatial 
mode, leading to efficient fiber coupling. The source is 
suitable for schemes that combine polarization and time 
coding. We have shown how distinguishability between 
photon-pairs is introduced for this type of colinear source, 
due to a special kind of chromatic two-photon dispersion. 
We have derived and analyzed the output state of SPDC 
for this case, with the goal to cancel the decoherence and 
regain a pure state using an extra piece of birefringent 
crystal. We have determined the quality of entanglement 
for the reported setup using various measures, including 
the method of quantum state tomography, and we draw 
the conclusion that this is one of the brightest sources 
available for polarization entanglement in terms of Bell- 
inequality violation and production rates. 



where |"0oo) is the number state at time to and H(t) is 
the interaction Hamiltonian, 

L/2 oc oo 

H(t) = J dz J d V Jdxx {2) El +) E^El~ ] +H.c, 

-L/2 — oo -oo 

(A.2) 

displayed in a Cartesian coordinate system, r = xe x + 
ye y + ze z . There arc three interacting fields in the crys- 
tal's volume, ignoring all higher-order terms (n > 3) of 
the non-linearity x ■ All three fields have the same 
polarization (ZZZ): 

E p +) = E e -*( fc o P Sp-r- Wp t+0 p ) (A. 3) 

4<-> - yd0 s yd Ws ^ s (c s )^ e J ^ s ="-^ t+ ^)at (WsiSs) 

(A.4) 

Jdfr yd^Ai^O^e^- 8 '—-'^^^^,^), 

(A.5) 
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where the pump field is classical and monochromatic so 
that we can replace Ep + ^ by E p + \ The plus sign denotes 
conjugation, i.e annihilation (+) or creation (-) of the 
state. We have also introduced the notation k = ks, 
where s = pe x + qe y + me z , is the unit length vector of 
k with components in each of the three dimensions |33j| , 
as defined by the coordinate system in Fig.^ The pump 
field is a plane wave propagating in the z-direction, s p = 
e z . For signal and idler, we sum over both frequency and 
angular modes, where a(uj, s) is the field operator, and 
A(lu) is the frequency amplitude of a Gaussian-shaped 
detector filter having the bandwidth AA (FWHM) and 
center wavelength A c (all wavelengths in vacuum). Via 
the relation to = 27rcn A /A, its form is given by 



APPENDIX: THE TWO-PHOTON FREQUENCY 
AND POLARIZATION QUANTUM STATE 

In this Appendix, we derive the quantum state of a 
single crystal in terms of frequency and polarization de- 
crees of freedom, using the interaction picture of SPDC 



The evolution of the number state vector is given by 

to+T 



exp 



1 

ih 



AtH{t) 



T 

to+T 



l^oo) 



1 + — / dtH(t) | |V>oo; 



(Al) 



(A.6) 



Each signal and idler photon is created with a random 
phase, 4> s and 4>i, respectively, which we need to sum over. 
The phase of the pump, p , is constant and arbitrary. 

For periodically poled materials, the spatial variation 
of the nonlinear index x' 2 ' has sharp boundaries, but we 
will simplify and make a sinusoidal approximation using 
the first term of an Fourier-series expansion of x' 2 ' : 



A 2 ) 



\ r — imK -r 



X2 2^ 



X2 /ie 



-iK- 



(A-7) 



where K = Ke z = 2n/Ae z and A is the grating period. 
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The Hamiltonian now takes the form 
H{t) = x 2 fiEa fdfafdfa Jdu s Jdwi 

x A s (uj s )Ai(uJi) 

x al(m s , s s )a{ (uj, Sj) 

s, Si 

L/2 oo oo 



x / dz J dy J dx 

—L/2 — oo — oo 



-i[Ak-(xe x +ye y +zez) — (uj s +uJi — u> p )t + s +0i — </> p ] 



H.c. 



where the mismatch vector is 



Ak — k s s s -\- k[S[ kopSp 



K 

Ak.e,. 



(A.8) 



(A.9) 



Following Eq. (|A.1(I . we now let the Hamiltonian 
undergo time evolution. The mismatch vector is also di- 
vided up into its x, y, and z components using Eq. 1|A.9J) . 
Hence, 



1 

ih 



dtH(t) 



X 2 .fi E o JdcJ B Jduji A s (uj s )Ai(oJi) 



s, S; 
L/2 oo oo 



Jdz Jdy Jdxt ' v ' \',.„..V.. ' 

-L/2 -oo -oo 

2tt T 

^JJ dfadfa J dt e -il(»s+"i-"v)t + 



H.c. 



<t>s+<t>i-4> P ] 
(A.10) 



The integration over the interaction volume, dx, dy, 
and dz, can now be easily carried out. There are three 
spatial integrals, of which two are the Fourier transforms 
of unity (dx and dy) and one is the transform of a box 
function (dz). The transforms turn into two S functions 
and a sine function, respectively. The time integral also 
turns into a S function of the three frequencies lo s , w;, 
and u> p . This is because we assume a monochromatic 
pump beam with infinite coherence length, which effec- 
tively leads to an infinite interaction time, T — > oo, even 
for short crystals. Motivated by the rotational symmetry 
of the emitted modes, we also change to a spherical coor- 
dinate system (see Fig. ^| , by replacing the summation 
over s with integrals over 9 Sl 9i, ip s and tpi. Furthermore, 
the only non-zero solution for the integration over the 
random phases, cp s and (pi, is for the phases to add up 
to a constant, yielding the relation (p s + (pi = <p v + C. If 



we let C = for simplicity, and drop some constants 
resulting from the integrations, we are led to 
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ih 



dtH(t) 



X2 h E o duj s duji A s (uj s )Ai(u}i) 
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n/2 

x J sin( 




k/2 2tt 2k 

d9 s J sin 9i dOi j d(p s J dtpi 
o oo 



x al(u) s ,0 s ,(p s ) a{(ui,6i,tpi) S(uj s + uji - lo p ) 



x S(Ak x ) S(Aky) L sine 
- H.c. 



-Ah 



(A.11) 



At this stage, we observe that fc s and ki each depend 
on to s and u>i, respectively. Motivated by the S- function 
in Eq. (|A.11|I . we let uj s = ujq s + e and u>i = — e, and 
make a series expansion of the k- vectors: 



h « fcoi - e 



dfcos 
duj 0s 

dfcpj 
do> i 



ko s 



>-0s 



(A.12a) 



1 



fcoi - = fcoi - e— • (A.12b) 



In a spherical coordinate system, we have p = sin 9 cos tp, 
q = sin 9 sin ip, and m = cos 9, and so the phase-mismatch 
vector components become 

Ak x = k s sin 9 S cos tp a + ki sin 9i cos tpi ps 0, 
Ak y = k s sin 9 S sin (p s + ki sin 9i sin pi ps 0, 



A/c z = k s cos 6* s + h cos 6*i 



«0p 



A' 



(A.13) 



where we have done a first-order approximation of sin# 
and cos# for small angles, meaning that we consider only 
plane waves, and where the last component is simplified 
using the phase-matching condition for the forward direc- 
tion, fc 0s + fcoi ~ ^0p + A = 0, together with Eq. (|A.12(1 . 
Thanks to Eq. (|A.13fl . we can now trivially perform the 
integration over the spatial modes d9 s , d9i and dtp, which 
finally leads to the following compact expression 



1 



dtH(t) 



1 



:X 2 fl E 
in 

x JdeA s (e)Ai(e)at(e)al(e) 



x L sine 
- H.c. 

de 17(e) at(e)a}(e) 



H.c. 



(A.14) 



In summary, Eq. (| A. 1 f> . via Eq. I|A.14(I . has helped us 
find the frequency and polarization state generated in 
one crystal, which we will write in the form 



^JdeU(e) 



\Xzz) 



(A.15) 
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where U(e) is defined by Eq. (|A.14|) . and where is a normalization constant, such that |-g JdeU(e)\ 2 = 1. 

Here, e represents the frequency mode and Xzz re P re ~ 
(Jde\A s (e)Ai(e)\ 2 smc 2 [Le(n^ s - n^ i )/2c]) 1 ^ 2 sents the polarization mode along the Z-axis. 



B = 



h( X2 fiE L)- 1 

(A.16) 
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